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Abstract 

Although many random-phase approximation (RPA) calculations of the Gamow-Teller (GT) 
response exist, this is not the case for calculations going beyond the mean- field approximation. We 






> . apply a consistent model that includes the couplmg of the GT resonance to low-lying v.brat.ons, 
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to nuclei of the fp shell. Among other motivations, our goal is to see if the particle-vibration 
coupling can redistribute the low-lying GT^ strength that is relevant for electron-capture processes 



^D ' in core-collapse supernova. We conclude that the lowering and fragmentation of that strength are 



consistent with the experimental findings and validate our model. However, the particle-vibration 
coupling cannot account for the quenching of the total value of the low-lying strength. 
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I. INTRODUCTION 

The Gamow- Teller (GT) resonances are among the clear manifestations of nuclear collec- 
tive motion. In (p,n) or (^He, t) reactions studied at small angles or zero angular momentum 
transfer, when the target nucleus is perturbed by an external field proportional to 

OGT-=E^C0-r-«, (1) 

i=l 

a broad peak systematically appears in all the nuclei that have been studied in the last 
three decades [l[. The position and the magnitude of such a peak can be explained only by 
assuming a coherent effect among particle-hole (p-h) excitations in the spin-isospin channel. 
Accordingly, much effort has been spent in trying to understand if the systematic study 
of the GT response can well constrain the nuclear interaction in the spin-isospin channel [2- 
|7|, which will be discussed in this work as well. From the nuclear structure point of view, 
another interesting phenomenon is the so-called Gamow- Teller quenching. The difference 
between the total strength associated, respectively, with the Oqt- and Ogt+ operators 
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[the latter being of course analogous to the former one but proportional to t+(z)], has t 
model-independent value given by 3(A^ — Z), which is the well-known Ikeda sum rule 
and will be denoted below by mo(GT~) — mo(GT+). Only about 60%- 70% of this value has 
been found around the main GT peak [9j]. In ^°Zr, it has been possible to assess the relative 



importance of different possible quenching mechanisms 10|]. In principle, the quenching may 



be due either to coupling with configurations that are more complicated than the simple p-h 
ones [ul, ll2| , or to coupling with internal degrees of freedom of the nucleon as with the A 
particle. Since 93% of the Ikeda sum rule is found below 50 MeV jlO], not much room is left 
for the latter mechanism. 

At the same time the Gamow- Teller and the other spin-isospin charge-exchange nuclear 
transitions do play a relevant role for particle physics [l3|, UJ^ and astrophysics 15|. We 
will not discuss much their relevance for single and double /3 decay, but we will focus on 
the connection between GT transitions and the evolution of massive stars at the end of the 
last hydrostatic burning phase. In fact, if the mass of the iron core exceeds the so-called 
Chandrasekar mass, then the pressure of the degenerate electron gas is not sufficient to make 
the system stable against gravitational collapse. In this scenario, one of the main processes 
that governs the subsequent evolution of the star until an eventual supernova explosion is 



the electron capture both by free protons and by core nuclei in the iron region [I5l4l8|. 

Extensive tables of electron capture rates, based on realistic estimates of electron capture 
cross sections, would be needed for supernova simulations. At the electron energies of inter- 
est, typically less than 30 MeV, electron capture is dominated by the GT+ contribution in 
nuclei around the mass region of Fe 
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-22l|. When an accurate experimental determination 
of the corresponding strength is missing, one is obliged to resort to theoretical frameworks 
such as the shell-model (SM) and the mean-field (MF) or energy-density functional (EDF) 
based schemes. EDF-based calculations of the electron capture process are worthwhile for 
several reasons. They are less demanding than SM calculations and are not limited to sd or 
fp shells. A cross-comparison with the SM can help in reducing the uncertainties associated 
with the implementation of the EDF-based scheme, and, in a complementary way, the same 
comparison is certainly instrumental to validate different EDF parameter sets, in the path 
toward a universal EDF. 



In Ref. [21I], the electron capture cross sections on several nuclei have been calculated 
by using a self-consistent implementation of the Skyrme Hartree-Fock plus random-phase 
approximation (HF plus RPA) model. This model, well known and widely used for many 
years, has been for the first time extended to finite temperature. The results have been 
compared with SM results from Refs. 
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The two models predict cross sections that are not very different at energies around 
10 — 20 MeV. At high energies, the SM predictions lie below the HF-RPA ones, whereas at 
low energies (below about 10 MeV) the situation is reversed. In fact, the main drawback 
of HF-RPA seems to be the fact that it predicts a too-high threshold energy for the GT^. 
Another feature of the finite temperature HF-RPA calculations is a certain spread of the 
results associated with different parameter sets. 

The present work is aimed to increase understanding and try to overcome these open 
problems. In the simpler zero-temperature case, where one can compare with experimental 
measurements of the GT+ strength, we first study the underlying reasons for the sensitivity 
of this quantity to the choice of the Skyrme parameter set. Then, we analyze the problem 
of the threshold energy and the fragmentation of the GT strength. The Skyrme functionals 
are well known to be unable to describe the single-particle states around the Fermi energy: 
they tend to predict a too-large single-particle gap, and cannot by definition reproduce 
the fragmentation of the single-particle strength. As the particle- vibration coupling (PVC) 



approach is one possible improvement in this respect, we will study to which extent a model 
based on RPA plus particle-vibration coupling can better reproduce the GT strength in the 
nuclei of interest. 

The plan of our paper is therefore the following. We provide first a description of the 
basic formalism in Sec. [ITlas well as of the numerical input in Sec. IIIII The results concerning 
the sensitivity to the choice of the Skyrme set are discussed in Sec. IIVI and the results of 
RPA plus particle-vibration coupling are discussed in Sec. |Vl We mainly concentrate on the 
nucleus ^°Ni, as a typical system in the mass region of interest. Meanwhile, we also show 
some results for other nuclei such as ^^Ni, which has been object of a recent experimental 



study |23|. We present conclusions in Sec. |VT1 



II. SKETCH OF THE THEORETICAL FRAMEWORK 



In the present work, we use the same basic formalism already employed in Ref . [2J] . The 
main difference is that the present PVC calculation is more consistent as we discuss below. 
The first step is a HF plus RPA calculation of the Gamow- Teller strength. This is done in 
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the same way as in Ref. J5J. The second step consists of implementing the coupling with 
vibrations. 

We start from the solution of the HF equations for a given nucleus {A, Z) using a Skyrme 
two-body interaction. Then, we set up a discrete basis of both proton-neutron and neutron- 
proton p-h configurations. The continuum is discretized by putting the system in a box and 
requiring vanishing boundary conditions for the wave functions at the surface of this box. 
We write the RPA matrix equations and solve them in the model space of the p-h discrete 



configurations, which is called Qi. These equations are well known from textbooks 25 1 
and we shall not discuss them further here. The RPA provides, as a rule, an accurate 
description of the centroid of giant resonances and the fraction of the energy weighted sum 
rule exhausted by the mode. 

Giant resonances have a quite large damping width. In order to explain it, one needs 



to go beyond the RPA. As discussed, e.g., in Ref. 2J], the width of giant resonances is 
mainly made up of an escape width (P^) and a spreading width (P"^). These are due to 
different mechanisms, namely to the decay through particle emission and to the coupling 
with more complicated states of the nuclear spectrum, respectively. The RPA can reproduce 



the escape width if proper couphng to the continuum is implemented, but cannot account 
for the spreading width. 



Here we will follow the formalism of Ref. [2J], and only focus on the spreading width 



which is the leading damping mechanism. The processes in which the energy and angular 
momentum associated with the vibrational nuclear motion are distributed among more com- 
plex internal degrees of freedom do contribute to the spreading width. In order to account 
for these effects we need a subspace, which we shall call Q2 and which is built with a set 
of "doorway states." We denote these doorway states by |A^), and make a physical choice 
of them in terms of states made up with a p-h excitation coupled to a collective vibration. 
Note that we use the same Skyrme interaction to calculate both the GT states and all the 
ingredients to build up the space Q21 namely we employ Skyrme-HF single-particle states 
and vibrations that are calculated consistently in the RPA with the same Skyrme set. In 
this respect, our scheme is parameter free. 

We use the projection formalism to restrict our effective Hamiltonian to the subspace Qi 
and make the calculations feasible. After truncation of higher-order couplings, the effective 
Hamiltonian reads 

n{uj) = Q,HQ, + W^iu) = Q,HQ, + Q.HQ^ ^ l^ , . Q2HQ1, (2) 

uj - Q2tiQ2 + ie 

where u is the excitation energy. This energy-dependent, complex Hamiltonian has complex 
eigenvalues whose imaginary parts originate from the coupling to the more complicated 
configurations. We shall denote these complex states as \v) in what follows. 

In practice, we first diagonalize the RPA Hamiltonian QiHQi and obtain the complete 
basis made up with the RPA states. The creation operators Ol of the states |z/) can be 
expressed as a linear combination of the RPA creation operators, namely 

Ol=Y. ^i'^Oi - Fi^^Ol (3) 

a;„>0 

where 0\ and O]^ are creation operators of the RPA states \n) lying at positive energy cu^, 
and of the states \n) at the corresponding negative energy — u;„, respectively. Then the 
eigenvalue equation for the effective Hamiltonian ([2]), that is, 

[n,ol] = {n^-i^)ol (4) 




can be cast in matrix form on the RPA basis as 

V + A,{UJ) AM W^ , .^ .^.^ - , ,,, 

-AM -V-Ai{uj) 

In this latter equation "D is a diagonal matrix with the positive RPA eigenvalues, and the 
Ai matrices contain the elements associated with the second term of Eq. ([2]) denoted by 

lv + Ai{uj) A(w) ^ 
W^. The matrix is complex and symmetric, as it can be seen 

\^ AM V + Ai{u)] 

from their explicit form provided below. The orthogonality and normalization relations of 
the eigenvectors are 

^^H^(.')_ir„HF„(^')^5^^,. (6) 

n 

Our goal is to extract from the diagonalization of the effective Hamiltonian ([2]) a physical 
observable. In particular, we are naturally interested in the response function associated with 
an external operator such as the GT operator defined in the Introduction. This response 
function can be written as 

R{uj) = (0|6J.T ^7^TX-^gt|0). (7) 

u —n[u) +irj 

We use a notation that includes both the GT~ and GT+ cases (one should note that O^rp. 
is equal to Ogt+j and vice versa). The strength function is related to Eq. ([7]) by the 
well-known relation 

5M = --Imi?(u;) = -hu^Y.{Q\doTW^-^-—r:. (8) 

where the squared matrix element of the transition operator appears, and not its squared 
modulus, due to the properties of the eigenvectors \v) which form a biorthogonal basis. 

We now provide some more details on the calculation of the matrix elements of W^. 
On the basis of the p-h configurations in Qi space on which the RPA has been solved, the 
matrix element of W^{uj) will be denoted as Wpj^^piy{(jj). If these matrix elements are known, 
then the matrix elements A on the RPA basis are obtained through a straightforward basis 
transformation, 

ph,p'h' 

iA,Un= E <y.'M^^r^>;^':u<,V.'(-^)C^^i'2'' (10) 

ph,p'h' 



ph,p'h' 

iA4,)mn = 2^ Wph,p'h'i.^Wph Yp'h' + ^ph,p'h'i~^)^ph ^p'h'- (12) 

ph,p'h' 

W'^ is clearly given by 

TV CO — UJn 

The matrix elements at the numerator can be evaluated using Wick's theorem. Wpf^p,f^,{u) 
turns out to be the sum of the four terms whose diagrammatic representation is shown in 
Fig. [T] and whose analytic expression is 

1 {p\\V\\p",nL){p'\\V\\f,nL) 



p",nL "^ y^'"' ' ^P" "^f^l ' "'I Jp 

1 {h\\V\\h",nL){h'\\V\\h",nL) 



u - {un + e.p'i -eh) + ir] j"} 



W\2) = 5,,,5,,„, E 



h^L ^-{^n- eh" + Cp) + IV 3i 

"^^u-iun + ep- eh') + 17] y j^, j^, L ' 
W\4) = E , , }~^T''^^'T\ ^ ,„ ^' ^' ^ \{p\\V\\p',nL){h\\V\\h',nL). 



(14) 



In the above formulas, p and h label always particle and hole states, respectively. The 
corresponding angular momentum and single-particle energy are given by ji and e^. jf is a 
shorthand notation for 2jj + 1. The phonon states are labeled by their angular momentum 
L (only natural parity states are included in our calculations) and by an additional index 
n. The (small) parameter rj is introduced to mimic couplings beyond the doorway-states 
approximation: this parameter is set at the value of 500 keV. 



Up to this point, the scheme is exactly the same as in Ref. 24[ . In Ref. 24j , the interaction 
at the particle-vibration coupling vertex \V in Eq. ( TT4l) ] was approximated by retaining only 
the Vo part in the following momentum-independent terms of the Skyrme p-h force (to and 
^3 terms), 

V,, = V^^{r, - r^) + V^'^dir, - r^)^, ■ ct^, 

V,^ = Vf5{r^ - ra) + Vf5{r^ - r<,)cT^ ■ cr^. (15) 
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FIG. 1: Diagrammatic representation of the four terms whose sum gives the matrix element 
W\ ,f^,. The analytic expressions are shown in Eq. ([1^ 



The functions Vq and V^- depend only on the radial coordinate r, and their detailed expres- 
sions are 

Vo'%r) = ho{l - xo) + ^hia + 2)(a + l)p"(r) - ^t,{x, + ^)p-{r) 

+^ha{l -«)(! + 2x3)p"-'(r)p24r) - ^^3(2x3 + l)ap"-i(r)p_(r), 



Vo'%r) = ito(2 + xo) + -^t3(a + 2)(« + l)p"(r) + -^t3(x3 
Z lb iz 

+lt3«(l - a)(l + 2x3)p"-2(r)p2_(r), 
K^'^(r) = ho{xo-l)-^t,{l-Xs)p-{r), 
Vf{r) = Itoxo + ^hx.p^ir). 



:)P°W 



(16) 



In these formulas, q and q label different charge states of two p-h pairs, where q ^ q. 
P- = Pn — Pp stands for the density difference of neutron and proton, while p = Pn + pp is 
the total nucleon density. The other symbols are the standard Skyrme parameters. In the 
present work, we include the whole central Skyrme p-h force, i.e., not only all the to and ^3 
terms, but also the ti and ^2 terms with the following form: 

1 



V = -ti(l + xiP.)[P^5(ri-r2) + 5(ri-r2)P'](l-P.P.) 

+t2(l + X2Pa)P' ■ 5{ri - r2)P(l + P.Pr), 



(17) 



where P = ^(Vi — V2), and P' is the Hermitian conjugate of P (acting on the left). 
Pa = 1(1 + CTi ■ cr2) and Pr = ^(1 + ti • r2) are the spin-exchange and isospin-exchange 



operators. This is consistent with the findings of Ref. 26|, where it 



las been shown that 



the momentum-dependent terms have a non-neghgible effect. In Ref. 26| the expression of 
the reduced matrix elements of V is derived. We report the result here, that is, 



{t\\V\\j,nL) = V2LTTY: [x;tVShjp) + (-)^+^'^-^''>;f ^l(^PJ/^) 

ph 

where the RPA amplitudes appear and Vl is the p-h coupled matrix element. 



(18) 



Viiihjp) = J2 



\3j-^j+3h-m.h 



all rn 



{jifnijj -nijlLM) {jpmpjh-mh\LM) {jiirii, jhirih \V\jjmj , jpUip) . 

(19) 



We briefly recall, for the reader's convenience, the relationship between the present ap- 
proach in which the full central part of the Skyrme interaction is considered in the PVC 
vertex, and the approximate one that was used in the past. If, as in Ref. 2J], only the 
Vo(r)5(ri — r2) part in to and ts terms is included in the interaction V, the p-h coupled 
matrix elements have the simple form 






— {^\\YL\\J){p\\YL\\h)j -K,'' (r)w.(r)«,(r)«p(r)«,(r). 



(20) 



where the radial part of the s.p. (single particle) wave functions, written as ipnijm 
^nh i^) j-y^ (g, ■\^^^2]jm, has been introduced. Thus, 



MlvunL) = V2LTiY.{x;k +Vk)ydih3p) 



ph 






nL ' 



(21) 



where q labels the charge of the states i and j, and the neutron or proton radial transition 
density of the state \nL) has been introduced. This is 



spi!kr 



V2L + ^pt^; 



j:ix;H' + Ypf){p\\yL\\h) 



Up{r)uh{r) 






(22) 



According to the further approximation used in Ref. 2J], i.e., namely that a collective 
state should have mainly isoscalar (or isovector) character so that its isovector (isoscalar) 
transition density can be neglected, Eq. (12T|) takes the even simpler form 



{i\\V\\3,nL) = {t\\YL\\j) / drV,'{r)u,ir)u,{r)5pi 



■iLi 



(23) 



which is exactly the form adopted in Ref. 24 1. 



Finally, we mention that, exactly in the same way as in Ref. 2J|, we introduce an isospin 
correction in the matrix elements of W'^ in the T_ channel. In fact, it must be noted that 
in this case the coupling with the intermediate states |A^) can mix states with different 
isospins. These intermediate states, which contain a phonon plus a proton particle and a 
neutron hole, do not have pure isospin. They can be written as 

|iV) = c_i|iV;To-l,To-l)+Co|iV;To,To-l)+c+i|iV;To + l,To-l), (24) 

where Tq = {N — Z)/2, and the coefficients q are simply Clebsch-Gordan coefficients, 

c_i = (1 - IToTolTo - l,To - 1) = (2To - l)'/V(2To + l)^/^ 

Co = (l-lToTo|To,To-l) = -(To + l)-i/2^ 
c+i = {l-lToTo\To + l,To-l) = l/[{To + iy/\2To + iy/% (25) 

(this writing is justified by the fact that our phonons, as mentioned below, are to a good 
approximation isoscalar phonons). The GT~ resonance has isospin quantum numbers 
\T,Tz) = \Tq — 1, To — 1), and its coupling with states of different isospin should be for- 
bidden by the nuclear Hamiltonian. We impose that it is strictly forbidden (this amounts 
to neglecting Coulomb effects in the residual interaction) and, therefore, we project out the 
isospin component with the same value of the GT~ resonance in the intermediate states, 
i.e., the component T = To — 1. Correspondingly, in the W'^ matrix elements, an isospin 
correction factor |c_ip = (2To — l)/(2To + 1) is added. 

III. NUMERICAL DETAILS 

In the present work, the HF equations are solved in coordinate space on a radial mesh 
whose size is 0.1 fm, within a box of 21 fm. The p-h configuration space used for the GT 
RPA calculation includes all hole states and particle states (discretized in the mentioned 
box) up to an upper cutoff -Ecut = 100 MeV. The results are fully converged in this way 
(for example, the Ikeda sum rule of ^°Ni with the interaction SGII 27| reaches 11.998 at the 
RPA level with this cutoff). 

To build the model space of the p-h pairs plus phonon doorway states, needed for W^, 
the energies and reduced transition probabilities of the most collective phonon modes with 
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TABLE I: Properties of low-lying phonons in Ni. The theoretical results are calculated with the 



Skyrme interactions SGII 
Ref. 1301. 



m SLy5 



28l | , and SkM* [29]. The experimental data are taken from 



State Theory Experiment 

Energy B{EL, 0^ L) Energy B{EL, -> L) 

[MeV] [e^ fm2^] [MeV] [e^ fm^^] 

SGII SLy5 SkM* SGII SLy5 SkM* 

2+ 1.820 1.157 2.212 3.611 x 10^ 3.562 x 10^ 2.457 x 10^ 1.333 8.780 x 10^ 

3^; 4.897 5.727 4.901 1.576 x lO'' 1.643 x 10^ 1.150 x lO'' 4.040 

if 1.902 2.343 3.469 1.007 x 10^ 7.257 x 10^ 6.365 x lO'' 2.506 

spin and parity 2+,3^,4+ are calculated with the same energy cutoff. With this choice in, 
e.g., ^°Ni, the energy- weighted sum rules (EWSRs) satisfy the double commutator values by 
about 98%. Non-charge-exchange RPA is implemented in exactly the same way as charge- 
exchange RPA, that is, with the same numerical input. 

To minimize violations of the Pauli principle, and be consistent with the very idea of 
particle-vibration coupling, only phonons which absorb a fraction of the total isoscalar or 
isovector strength larger than 5% (and with energy less than 20 MeV) are included in the 
model space Q2- The properties of the low-lying 2^^, 3~[, and 4^^ states in ^°Ni calculated with 
three different Skyrme interactions, namely SGII, SLy5, and SkM*, are shown in Table HI 

We finally point out that we have checked, in the present calculation of Gamow- Teller 
resonances, that the Ikeda sum rule is still satisfied at the level of RPA plus particle- vibration 
coupling. For instance, in the case of ^°Ni calculated with the Skyrme interactions SGII, 
SLy5, and SkM*, mo(GT") - mo(GT+) = 11.87,11.89, and 11.83, respectively, when the 
strength is integrated up to 45 MeVfor ^^Ni, the sum rule value is -0.10 for the interaction 
SGII, -0.07 for SLy5, and -0.11 for SkM* when the strength is integrated up to 45 MeV; and 
for ^°*Pb it is 129.29 for the interaction SGII when the strength is integrated up to 70 MeV. 



11 



TABLE II: Correspondence between numbers and interactions in Fig. [21 Fig. [31 and Fig. [H 

Number 12 3 4 5 6 7 8 9 10 11 12 13 

Interaction SLy5 SLy4 SkP LNS SkI3 SkM* Sk255 SIII BSkl7 SGII SkO' SkI4 SkO 

Reference [28] [28][31j[32][33][2^ [3^[35j[36l [2;a[3r][33j[3r] 



IV. SENSITIVITY OF THE RPA RESULTS TO THE CHOICE OF THE SKYRME 
SET 

We have first tried to understand the sensitivity of the GT energies to the choice of the 
Skyrnie interaction used for RPA, and the underlying reasons. 

The main components of the GT response are the transitions between spin-orbit partners. 
In a nucleus such as ^''Ni, which is our benchmark in this study, the GT^ (GT+) response 
is expected to be dominated by the single transition Z//7/2 — ?■ T^f^i2 {j^fii2 -^ ^f^i2)- In such 
a case, we can approximately write the energy of the GTR as 

^RPA = ^unper + ^res ~ AE,, + V. (26) 

The first term -Eunper is the unperturbed p-h energy, and can be associated with the spin- 
orbit splitting l^Eis whereas the second term ii^res is the energy shift induced by the residual 
interaction and can be written in terms of its matrix element V . The matrix elements of 
the residual interaction in the spin-isospin channel are as a rule repulsive, so the two terms 
of Eq. (!26l) have the same sign. 

In the case of a Skyrme calculation, the first term is controlled by the spin-orbit parameter 
PVo (we leave aside the nonstandard Skyrme sets that have two parameters in the spin-orbit 
part of the energy functional, and those which include the so-called J^ terms). The second 
term depends instead on the strength of the residual interaction in the spin-isospin channel. 
In general one imagines that it is mainly associated with the well-known Landau parameter 
(7q, but we should not forget that g'^ plays a role as well (these are the only two Landau 
parameters that do not vanish for a zero-range interaction). In Ref. p], by studying the 
strength of the GT~ resonance in several nuclei including Sn isotopes and ^°^Pb, we have 
concluded that this strength is sensitive to (^g and f^J, and that values around 0.45 and 0.5 
are preferable for these parameters. 

In the present study, we focus first on the energy position of both the GT~ and GT"*" 
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FIG. 2: (Color online) The Gamow- Teller energy calculated in the RPA (E^rpa), the unperturbed 
energy of the main p-h components (-Eunper, see the text), and the associated shift ii^res = -E'rpa — 
-E'unper are displayed as functions of the Landau parameter g'^ for nucleus ^"^Ni. The GT^ and GT"*" 
cases are in panel (a) and panel (b), respectively, and the correspondence between numbers and 
Skyrme sets can be found in Table [III 

peaks. In Fig. El the following quantities are displayed, as a function of (^g: the main peak 
of the Gamow- Teller resonance calculated in RPA, the unperturbed Gamow- Teller energy 
of the single p-h configuration {7if5/2){^f7/2)~^ for GT~ or {v f 5/2) {t^ 17/2)^^ for GT+, and 
the shift -Eres = -S'Rpa — -E'unper- Surprlsiugly, the energy shift is fiat as a function of Qq. we 
have found that this is due to a strong cancellation between the contributions associated 
with q'q and (7^. The staggering of the GT~ main peak is associated with the staggering in 
the unperturbed energy (which is obviously uncorrelated with the Landau parameters), and 
it follows it quite closely. The same can be said of the GT+ peak; in this case the absolute 
value of the unperturbed energy is, however, smaller than in the previous case and smaller 
than the typical matrix element of the residual force. 

In Fig. |3] the RPA and unperturbed energies are again displayed as a function of the 
spin-orbit strength parameter Wq. Only a subset of all the studied interaction has been 
chosen for this figure. The dashed and dash-dotted lines are linear fits of the points that are 
displayed only as guides to the eye. A quantitative correlation between the parameter Wq 
and the GT energy can be neither inferred from the numerical results nor expected, since 
the unperturbed energy is associated with a neutron-proton transition. However, it is clear 
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FIG. 3: (Color online) The same energies -Erpa and -Eunper shown in Fig. [2] are here displayed as 
functions of Wq. The lines are linear fits that are, however, shown only for illustrative purposes. 
See the text for a detailed discussion. The experimental values for the energies (expressed with 
respect to the parent nucleus as in the theory) are indjcated by means of black solid lines and are 
taken from Ref. [38] for the GT~ case and from Ref. [39'] for the GT+ case. 

that there is a qualitative tendency of both unperturbed and RPA energies to grow as a 
function of Wo, and that these energies are more sensitive to the single-particle spectrum 
than to other parameters such as (7q or g[. 

In the same figure we also dis play, by means of horizontal lines, the experimental values 



of the energies taken from Refs. 



38 



Since the particle-vibration coupling produces, as 
a rule, a downward shift of the resonance energies (see the next section), we are inclined to 
prefer, for further considerations, the sets marked as 2, 6, 8, and (at least looking at the 
GT~ case) 10. Set 8 is SIII, which is a force that does not have very satisfactory properties 
as a whole (for instance, the associated value of the nuclear incompressibility is too large). 
Set 2 is SLy4, but SLy5 provides quite similar results. Sets 6 and 10 are, respectively, SkM* 
and SGII. We will discuss in the next section some results obtained by using SkM*, SLy5, 



and SGII. SkM* and SGII have reasonable values for the Landau parameters 



5|- 



In Ref. p] it has been found, however, that the parameters q'q and g[ are not irrelevant for 
the GT properties, in keeping with the fact that they are correlated with the GT" strength 
in 208p]^ g^g ^g^ g^g jj^ g]^ isotopes. We have checked here what happens for the GT~ energy 
in ^°^Pb, and our results are shown in Fig. |H Although the single-particle properties still 
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FIG. 4: (Color online) The same quantities already displayed in Fig. [2] are now displayed in the 
case of ^°^Pb. Here, the unperturbed energy -Eunpcraver is the weighted average of the two main 
configurations 1^113/2 — ^ '^*ii/2 ^''^'^ ^^11/2 ~^ ^^9/2- 

play a relevant role - as can be seen from the staggering that is similar to that in the figure 
for ^°Ni - here there is a clear tendency of the shift associated with the residual interaction: 
it is not flat, rather it grows to some extent with Qq. We have found that in fact there is 
less cancellation between the terms associated with g'^ and g[. Moreover, in the schematic 
model, the RPA energy is known to scale as 



E- 



RPA 



E, 



unpcraver 



nV, 



(27) 



where n is the number of p-h configurations, -Eunperaver is the weighted average of the unper- 
turbed energies of the n configurations contributing to the RPA state. This number is not 1 
here, as it was in the case of ^"^Ni. The role of the residual interaction must be larger here, 
accordingly. 

In conclusion, the GT energy is, within the Skyrme framework, very sensitive to the 
single-particle spectrum. This is especially true in nuclei of the fp shell (although we have 
discussed only ^°Ni in some detail, we have reached similar conclusions from the analysis of 
56,58]\]-j ^-j^^ ^"^'^^Fe). There is some kind of tendency for both the GT~ and GT+ energies 
to grow as a function of the spin-orbit parameter Wq. Our analysis has shown that a few 
Skyrme sets may be suitable to be employed for PVC calculations with the expectation that 
the experimental energy can be well reproduced. 
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FIG. 5: (Color online) Ganiow- Teller strength distributions calculated with the Skyrme interac- 
tion SkM* for the nucleus ^'^Ni. The blue dashed discrete lines denote the RPA strength (with 
dimensionless units), and the black solid lines represent the distribution (with units of MeV~^) 
calculated by the RPA+PVC model. The experimental data for the case of GT~ [38] [with units 
of mb/(sr MeV)] and GT+ [39|] (with units of MeV~^), displayed by either a red dash-dotted line 
or red points, are discussed in the text. 

V. RPA WITH PARTICLE VIBRATION COUPLING 



In the present section we discuss the results of the calculations that include the coupling 
with vibrations. We use the interactions SkM*, SLy5, and SGII, following the discussion in 
the previous section. The goal is to see if, and to what extent, the coupling with vibrations 
decreases the energies (in particular the GT+ threshold energy) and changes the distribution 
of the strength. Of course we are also interested in seein g if RPA plus P VC can reproduce the 
total resonance width. In our previous work for ^'^^Pb 2J|, we have seen that the coupling 
with vibrations can account for most of the spreading width of the GTR. 



A. Energies and widths of the GT states in ^''Ni 

The Ganiow- Teller strength distributions as a function of the energy (with respect to the 
parent nuclei) are displayed in Figs. [5l El and [7] in the case of the Skyrme interactions SkM*, 
Sly5, and SGII, respectively, for the nucleus ^°Ni. The blue dashed discrete lines denote the 
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RPA strength, and the black sohd hnes represent the strength distribution calculated by the 
RPA+PVC model. From all figures, it is clear that the inclusion of phonon coupling leads to 
two significant effects for both the GT~ and the GT"*" case: there is a downward shift of the 
excitation energy, and the development of a spreading width. Especially in the high-energy 
region, a single RPA GT^ peak spreads into a wide resonance covering an energy region of 
the order of 5 MeV (or more). 

The development of a realistic spreading width allows, in principle, a full comparison of 
the theoretical results with the experimental data. The spreading width is caused by the 
coupling with the phonons, through the inaaginary part of W^. In the case of GT~, a detailed 



strength distribution is not found in Ref . J38j , so only the measured zero-degree cross section 
is displayed in Figs. \5\ Eland [7] by means of a dash-dotted red line (cf. the left panels). All 
three interactions produce two peaks in the low-energy part, and one main high- lying GTR. 
These qualitative features are in agreement with the experimental findings. However, the 
energies of the low-lying peaks are underestimated, and the GTR is not as broad as in the 
experiment. One can say we can account for a relevant fraction of the measured total width, 
yet the high- lying experimental cross section looks extremely fiat. 



In the case of the GT^ the experimental data from Ref. 39|] are indicated by points. The 
total width of the GT"*" distribution is about 3 MeV. The theoretical result obtained with 
the force SkM* agrees very well with experiment, as far as both the peak energy and the 
width are concerned. SLy5 also reproduces well the width, whereas in SGII one finds two 
peaks but still a fragmentation of the strength over a reasonable energy interval. It must 
be noted that the experimental data are multiplied by a factor of 3 in the right panels of 
Figs. [5l [6] and [71 We discuss this issue in the next subsection. 

B. The quenching problem 

As a result of the phonon coupling, we can say we redistribute the Gamow- Teller strength 
in a region that is several MeV wide. However, this does not produce a real quenching. In 
the case of ^"Ni with the force SGII, up to the energy of 20 MeV (10 MeV), the integrated 
GT" (GT+) strength calculated by RPA+PVC is 19.4 (8.1). For comparison, the RPA 
value is 20.8 (9.3). Consequently, the strength is quenched by 7% (13%) compared with 
RPA but it is larger than the experimental data. In experimental results, one finds a total 



17 



strength equal to 7.2 ± 1.8 in the case of the GT [38| (up to 21 MeV in the parent nucleus) 



and 3.11 ± 0.08 in the case of the GT"*" 39|] (up to 10 MeV in the parent nucleus). In the 
calculation of RPA+PVC with interaction SGII, only 15% of the sum rule shifts to energies 
above 20 MeV. 

SkM* and SLy5 give similar (i.e., small) quenching effects as SGII. In the case of SLy5, 
one finds an integrated strength of 23.4 up to 20 MeV for the GT^ case after coupling with 
phonons (instead of 24.8 given by RPA). For the GT+ case, the integrated strength up to 10 
MeV in RPA+PVC is 12.1 (instead of 13.0 from RPA). 11% of the total sum rule is shifted 
to the energy region 20 — 45 MeV. In the case of SkM*, the integrated strength for GT~ 
(GT+) from RPA+PVC is 20.3 (9.5) up to 20 MeV (10 MeV), instead of 21.9 (10.4) from 
RPA. 15% percent of the sum rule is shifted to the energy region 20 — 45 MeV. 

There are probably two missing effects that can explain the lack of quenching in our 



calculations (let alone the small, or negligible, coupling with the A isobar). In Ref. [ll| a 
large shift of strength to high energies is found by coupling with uncorrelated 2p-2h doorway 
states. These extend up to higher energy, as compared with the Ip-lh plus phonon doorway 
states. Consequently, although less important for the spreading width, they can account for 
the shift of strength above 20 MeV. In that work, the tensor force has been found to play a 
role as well. Even in the limited Ip-lh space, the tensor force has been found to be capable 

Q 

of shifting strength at high energy, although in smaller amounts [7|. 

C. Microscopic origin of the downward shift induced by phonon coupling 

In general, it can be expected that the coupling with phonons leads to a downward shift 
of the resonance peaks, and that the low-lying vibrations are the most effective. We analyze 
here in detail the effects in the case of the GT+ peak of Fig. [71 

The total shift in the peak energy for the GT+ channel is —2.1 MeV, and the 2+, 3~, 4+ 
phonons give partial contributions equal to —1.3, —0.8, and —0.1 MeV, respectively. For 
every multipole, the lowest phonon, which is the most collective one, plays the most im- 
portant role, giving 85% (-1.1 MeV) and 75% (-0.6 MeV) of the shift for the 2+ and 3" 
case, respectively. The diagonal matrix element W^ip/i, when the ph configuration is the 
dominant component of the GT+ peak [z^l/5/2(7rl/7/2)~^], is identified as the dominant term 
responsible for the energy shift. In the case of coupling with the first 2+ phonon, the first 
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[Fig. [T] (1)] and second [Fig. [T] (2)] diagrams of W^i^p^^ give similar negative contributions 
to the energy shift both with values of about —1.3 MeV, while the third [Fig. [1] (3)] and 
fourth [Fig. [T] (4)] diagrams cancel them by about 70%. These numbers are calculated at 
the value of u at which the GT peak lies if we couple the RPA state with the 2f state 
only, namely u = 1.4 MeV. In the first and second diagrams of Fig. [H the 2"*" phonons can 
only couple to the intermediate states with negative parity. It turns out that the particle 
states z^l/5/2 and i'2pi/2 are the most important intermediate states for the first graph, and 
provide contributions accounting for 62% (—0.8 MeV) and 23% (—0.3 MeV) of the total 
value of the first diagram (—1.3 MeV), while the intermediate hole state vrl/7/2 is the most 
important one for the second graph, giving almost all the contribution to the total value 
of the second diagram (—1.3 MeV). In the case of coupling with the first 3~ phonons, the 
first and second diagrams still give similar negative contributions to the energy shift (both 
about —0.3 MeV, these values being calculated at the value of u at which the GT peak lies 
if we couple the RPA state with the 3j~ state only, namely u = 1.9 MeV); however, the 
third and fourth diagrams are zero in the case of Wpj^pj^ when ph is z/l/5/2(vrl/7/2)~^, due to 
parity conservation. In this matrix element, the 3~ phonons only couple with intermediate 
states with positive parity. The states i'lgg/2, ^^97/21 ^^2^5/2 are the most important particle 
states and the associated contributions to the first diagram of Fig. [1] are —0.05, —0.05, and 
—0.03 MeV (the total value of the first diagram is —0.3 MeV). Similarly, in the second dia- 
gram the important hole states are '7r2si/2, vrl(i5/2, '7rl(i3/2, and the contributions are —0.13, 
—0.08, and —0.05 MeV (the total value of the second diagram is —0.3 MeV). 

D. Results for the nucleus ^^Ni 

Since the GT^ resonance in the nucleus ^^Ni was recently measured [23|, we also present 
results for this nucleus, calculated with the three interactions SkM*, SLy5, and SGII (cf. 
Figs. [8|9| andHni respectively). Similarly to ^°Ni, the coupling with phonons shifts the peaks 
downward and produces a spreading width, yet brings a quite limited quenching effect. The 
energy shift caused by the phonon coupling is about 1 — 2 MeV for all three interactions. 
The energy calculated by means of SGII and SLy5 is somewhat lower than the experimental 
data, while the energy from SkM* agrees very well with them. In this latter case (SkM*) the 
strength spreads out in the same energy region as the experimental strength distribution. In 
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FIG. 6: (Color online) The same as Fig. [5] in the case of the Skyrme interaction SLyS. 
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FIG. 7: (Color online) The same as Fig. [5] in the case of the Skyrme interaction SGII. 

experimental results, there is a double-peak structure and each peak has a width of about 
2 MeV. In our calculations, we also obtain a similar structure in the strength distribution 
for all three interactions. For SGII and SLyS, the first peak has a width of about 1 MeV 
and the second peak has a width of about 2 MeV. For SkM*, the width is about 1.5 MeV 
for the first peak and 1 MeV for the second one. The big difference between our results and 
the experimental data is the absolute value of the strength. As in ^"Ni, the coupling with 
phonons does not induce a significant quenching (see the above discussion). In the energy 
region 10 — 24 MeV, after the coupling with phonons, the integrated strength becomes 
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FIG. 8: (Color online) Ganiow- Teller strength distributions calculated with the Skyrme interac- 
tion SkM* for the nucleus ^^Ni. The blue dashed discrete line denotes the RPA strength (with 
dimensionless units), and the black solid line represents the distribution (with units of MeV~^ 



calculated by the RPA+PVC model. The experimental data (with units of MeV ^) from Ref. 
are indicated by red points. 
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FIG. 9: (Color online) The same as Fig. [8] in the case of the Skyrme interaction SLy5. 

10.3, 12.8, and 10.8 (instead of 11.0, 13.4, and 11.5 calculated in RPA), for the interactions 
SGII, SLy5, and SkM*, respectively. The same value is 3.5 ± 0.3 ±1.0 from experiment 23 |. 
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FIG. 10: (Color online) The same as Fig. [8] in the case of the Skyrnie interaction SGII. 
VI. CONCLUSION 

In the present work, we have analyzed the features of the Gamow- Teller strength distri- 
butions in nuclei of the fp shell. These nuclei are of interest for astrophysical applications, 
but few have been also objects of experimental measurements. Our purpose is to improve 
existing models so that they can reproduce the experimental findings, whenever available, 
and become reliable for calculations of the electron capture rates needed for astrophysical 
simulations; our ambition is to avoid resorting to free, ad hoc parameters. 

We have first performed RPA calculations with many different Skyrme parameter sets. 
We have seen that the energy of the GT^ and GT+ peak is quite sensitive to the single- 
particle properties and in particular, to a large extent, to the spin-orbit strength parameter. 
The energy shift caused by residual interaction is less relevant to the Landau parameter gl^ 
for these medium-mass nuclei than in the case of, e.g., ^°^Pb, the reason being the smaller 
number of p-h configurations. 

Coupling with phonons is relevant to producing a more realistic strength distribution, 
characterized by a spreading width. The interaction SkM* does reproduce the peak po- 
sition and the spreading width, in the cases of the GT"*" resonance in ^"Ni and the GT~ 
resonance in ^^Ni, after the phonon coupling is taken into account. This is important for 
the astrophysical applications since the GT"*" energy is associated with the threshold energy 
for electron capture. Other interactions display nonetheless, in the same nuclei, realistic 
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strength fragmentation after the phonon couphng is taken into account. 

We omit a full description of the GT strength quenching. While in shell-model calcu- 
lations this is often accounted for by means of a free parameter, we refrain from making 
this choice here. Our model probably can be improved by including coupling with doorway 
states extending to higher energy, and by taking care of the tensor force as well. 
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